Quasi-one-dimensional antiferromagnetism and multiferroicity in CuCrC>4 
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The bulk magnetic properties of the new quasi-one-dimensional Heisenberg antiferromagnet, 
CuCrC>4, were characterized by magnetic susceptibility, heat capacity, optical spectroscopy, EPR and 
dielectric capacitance measurements and density functional evaluations of the intra- and interchain 
spin exchange interactions. We found type-II multiferroicity below the Neel temperature of 8.2(5) K, 
arising from competing antiferromagnetic nearest-neighbor (J nn ) and next-nearest-neighbor (Jnnn) 
intra-chain spin exchange interactions. Experimental and theoretical results indicate that the ratio 
Jnn/Jnnn is close to 2, putting CuCrCU in the vicinity of the Majumdar-Ghosh point. 

PACS numbers: 75.30.Et,75.40.Cx,75.85.+t 
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Introduction 



Ferroelectricity driven by magnetic ordering in so- 
called type-II multiferroics exhibits a high potential for 
technological applications. Switching ferroelectric po- 
larization by a magnetic field or magnetization by an 
electric field offers unprecedented applications in mod- 
ern energy-effective electronic data storage technology* 1 * 2 - 
However, the link of magnetic order and ferroelectric- 
ity in type-II multiferroics still remains an intriguing 
question^ - — To elucidate this issue, lately much atten- 
tion has been focused on the magnetic and magneto- 
electric (ME) properties of quasi-one-dimensional (ID) 
antiferromagnetic (afm) quantum chain systems, which 
exhibit incommensurate cycloidal magnetic ordering^ 
Such systems lose inversion symmetry and appear to 
be suitable candidates for multiferroicity. Incommensu- 
rate spin-spiral magnetic ordering occurs in magnetic sys- 
tems consisting of ID chains when the intrachain nearest- 
neighbor (nn) and next-nearest-neighbor (nnn) spin ex- 
change interactions ( J nn and J nn n, respectively) are spin- 
frustrated, as found for compounds with CuX2 ribbon 
chains made up of CuX 4 plaqucttcs, where X is a suit- 
able anion, e.g. oxygen or a halide. Current examples in- 
clude LiCuV0 4 , NaCu 2 2 , CuCl 2 £"ii It is typical that 
the Cu-X-Cu superexchange J nn is ferromagnetic (fm), 
the Cu-X. . .X-Cu super-superexchange J nn n is afm and 
larger in magnitud e 12 i 15 A cycloidal spin-spiral along a 
ID chain induces a macroscopic electric polarization, P 
oc ey x (Si x Sj), where e.y is the vector linking the mo- 
ments residing on adjacent spins Si and Sji^— 

In an ongoing effort to identify new quantum spin chain 
systems which potentially exhibit spiral magnetic order 
and ferroelectric polarization, we recently focused our at- 
tention on compounds crystallizing with ribbon chains, 
mainly those belonging to the CrV04 structure-type. 
The aforementioned structure-type features MO2 ribbon 
chains where M is a magnetic id transition metal. Such 



compounds were recently shown to exhibit exotic mag- 
netic ground-states*i2r— Here, we report on the magnetic 
and ME properties of another member of this structure- 
type, CuCr04. Our density functional calculations in- 
dicate J nn to be about twice as strong as J nnn putting 
CuCr04 in the vicinity of the Majumdar-Ghosh point 
for which the ground state can by exactly solved* 2 ^ This 
feature makes CuCr04 uniquely exceptional since all of 
the CuX 2 ribbon chain systems investigated so far ex- 
hibit fm J nn and afm J nnn spin exchange, where J mm is 
considerably larger in magnitude than J nn & 12 i 15 ? 4 We 
demonstrate that CuCr04 exhibits long-range afm or- 
dering below ~ 8.2 K, which is accompanied by a ME 
anomaly due to possible spin-spiral ordering in the Cu0 2 
ribbon chains. 



I. CRYSTAL STRUCTURE 

CuCr04 crystallizes in the CrV04 structure- typ o 25 : 26 
(SG: Cmcm, No. 63) with Cu 2 + (d 9 , S =1/2) and Cr 6+ 
(d°) ions. In the crystal structure of CuCr04, the axially- 
clongated CuOe octahedra share edges to form chains 
running along the c-axis (Fig. [TJa)). These chains are 
interconnected by Cr04 tetrahedra such that each Cr04 
tetrahedron is linked to three CUO4 chains by corner- 
sharing (Fig. DIb)). The x 2 -y 2 magnetic orbital of each 
CuOg octahedron is contained in the CUO4 plaquctte 
with four short Cu-0 bondsi^i Thus, as far as the mag- 
netic properties are concerned, CuCr04 consists of cor- 
rugated Cu0 2 ribbon chains running along the c-axis 
(Fig. QTa)). At room temperature the Cu. . .Cu distance 
is 2.945(2) A and the Cu-O-Cu Z is 98.1(1)°. 



II. SPIN EXCHANGE INTERACTIONS 

To examine the magnetic properties of CuCr04, we 
consider the four spin exchange paths defined in Fig. [T] 
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obtain the values of the spin exchange parameters, J n 
J nnn , Jx, and J r 2 as summarized in Table [I] 




(-4J nn -4J nnr ,-8J(-4J 2 )/4 (+ , 'J nn -4J nnn -BJ 1 -4Jj)y4 (-4J nn -4J nnn +4J z )/4 
[0,0] [-34.43,19.98] [3.SS, 2.81] 



b 

J- 



(-4J nn -4J nnn +8J r 4J 2 )/4 (+4J nnr ,-8J r 4J J/4 
[-2.96,-2.05] [-32,00,-19.73] 



FIG. 1: (Color online) (a): The crystal structure of CuCrCU. 
The (blue) octahedra are the CuC>6 units while the (green) 
tetrahedra are the CrC>4 units. The interchain spin exchange 
pathways Ji and J2 are also indicated, (b): A section of the 
CUO2 ribbon chain highlighting the edge sharing CuCU pla- 
quettes, with the nn J nn and nnn J nnn spin exchange pathways 
labeled. 



FIG. 2: Five ordered spin states used to extract the values 
of Jnn, Jnnn, Ji , and J2 , where the Cu 2+ sites with different 
spins are denoted by filled and empty circles. For each ordered 
spin state, the expression for the total spin exchange energy 
per 4 FUs is given, and the two numbers in square bracket 
(from left to right) are the relative energies, in meV per 4 
FUs, obtained from the GGA+[/ calculations with U e s = 4 
and 6 eV, respectively. 



the two intra-chain exchanges J nn and J nnn as well as 
the inter-chain exchanges J\ and J^. To determine the 
values of J nn , Jnnn, J i and J2, we examine the relative 
energies of the five ordered spin states depicted in Fig. [2] 
in terms of the Heisenberg spin Hamiltonian, 

H = —'^^JijSiSj, (1) 

where Jy is the exchange parameter (i.e., J nn , </nnn, J\ 
and J2) for the interaction between the spin sites i and 
j. Then, by applying the energy expressions obtained for 
spin dimers with N unpaired spins per spin site (in the 
present case, N = l)r& the total spin exchange energies of 
the five ordered spin states, per four formula units (FUs), 
are given as summarized in Fig. [2] We determine the rel- 
ative energies of the five ordered spin states of CuCr04 
on the basis of density functional calculations with the 
Vienna ah initio simulation package, employing the pro- 
jected augmented-wave method^ - — the generalized gra- 
dient approximation (GGA) for the exchange and corre- 
lation functional^ with the plane-wave cut-off energy 
set to 400 cV, and a set of 64 k-points for the irreducible 
Brillouin zone. To account for the strong ele ctron cor- 
relation associated with the Cu 3d state, we performed 
GGA plus on-site repulsion (GGA+?7) calculations with 
U e g = 4 and 6 eV for Cui22 The relative energies of the 
five ordered spin states obtained from our GGA+U cal- 
culations are summarized in Fig. [2J Then, by mapping 
these relative energies onto the corresponding relative en- 
ergies from the total spin exchange energies ) 27 ' 34 " — we 
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-85.8(0.43) 
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-8.6(0.04) 
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-6.00(0.05) 
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+31.1(0.16) 


+9.8 


+22.3(0.19) 


+10.7 


+ 12 


dew 




-43.2 




-38.8 


-56/-60 



TABLE I: Spin exchange parameters J nn , Jnnn Ji and J2 (in 
K) of CuCr04 obtained from GGA+C7 calculations with U c s 
= 4 and 6 eV. The left column for each U a s contains the 
theoretical results, while the values in the right column are 
the scaled theoretical results such that J nn n equals the ex- 
perimental finding, -27 K. The rightmost column summarizes 
the experimentally found spin exchange values. The final row 
show the Curie- Weiss temperatures of the scaled GGA+J7 
spin exchange parameters, calculated using the mean field 
expression; #cw = | '^2 i ZiJiS(S + 1), where Zi is the num- 
ber of neighbor with which a single atom interacts with the 
spin exchange Ji, and the experimentally observed values (see 
below) . 

The intra-chain spin exchanges J nn and J nn n are both 
afm and constitute the two dominant spin exchanges in 
CuCr04. The inter-chain parameter J2, connecting Cu 
atoms related by a translation along a, is fm and, depend- 
ing on the onsite repulsion parameter U e g, its magnitude 
amounts to 15 to 20% of the intra-chain spin exchange 
J nn . Ji, which couples adjacent spin moments which are 
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related by a translation along [110], is afm and compara- 
tively small. Therefore, to a first approximation, CuCr04 
can be described as a quasi ID Heisenbcrg magnet with 
nn and nnn spin exchange interactions, both being afm. 
Since these ID chains are connected by weak inter-chain 
exchanges (Ji and J2), long range ordering will eventu- 
ally take place at low temperatures. 



A. Experimental 



atom Wyckoff site 



y 



Cu 4a 

Cr 4c 

01 8f 

02 8g 



0.09( 8) 

0.3700(3) 0.25 0.90( 8) 

0.2652(5) 0.0320(9) 0.80(12) 

0.2326(7) -0.0198(6) 0.25 0.80(12) 



TABLE II: Atomic positional parameters of CuCrCU (SG: 
Cmcm) as obtained from a profile refinement of the x- 
ray powder diffraction pattern, collected at room tempera- 
ture. The lattice parameters amount to a = 5.4388(5) A, 
b = 8.9723(8) A and c = 5.8904(6) A. 



A polycrystalline sample of CuCr04 was prepared by 
separately dissolving equimolar amounts of anhydrous 
Coppcr(II)acetate and Chromium(VI)oxide in distilled 
water, similar to the recipe given by Arsene et al£^. The 
two solutions were mixed and boiled to dryness. The re- 
sulting powder was heat treated in air at a temperature 
of 150°C for 2 days. The phase purity of the sample was 
checked by x-ray powder diffraction measurements using 
a STOE STADI-P difffactometer with monochromated 
Mo-ifai radiation. The powder pattern was analyzed 
using the Rictvcld profile refinement method employed 
within the Fullprof Suited No other reflections besides 
those of CuCrC>4 were observed. 

Powder reflectance spectra of CuCr04 were collected 
at room temperature using a modified CARY 17 spec- 
trophotometer, equipped with an integrating sphere. The 
spectrometer was operated in the single-beam mode us- 
ing BaSC>4 as reflectance (white) standard. CuCrC>4 
powder was mixed with BaSC>4 in a volumetric ratio 
CuCr0 4 :BaS0 4 ~1 : 5. 

Temperature dependent electron paramagnetic reso- 
nance (EPR) spectra of a ~ 5 mg polycrystalline sam- 
ple, contained within an EPR low-background suprasil© 
quartz tube, were collected using ~ 9.5 GHz microwave 
radiation (Brukcr ER040XK microwave X-band spec- 
trometer, Bruker BE25 magnet equipped with a BHf5 
field controller calibrated against Diphenylpicrylhydrazyl 
(DPPH)). 

The molar magnetic susceptibilities, x mo i, of a poly- 
crystalline sample weighting ~ 84 mg were measured 
with various fields between 2 K and 350 K using a 
SQUID magnetometer (MPMS, Quantum Design). The 
raw magnetization data were corrected for the magneti- 
zation of the sample container. 

The specific heats, C p , of a powder sample weighting 
~ 2.4 mg were determined as a function of the tempera- 
ture and magnetic field with a relaxation-type calorime- 
ter (PPMS, Quantum Design) for the temperature range 
0.4 K to 50 K and magnetic fields up to 9 T. 

The relative dielectric constant, e r , was measured at 
a constant frequency and excitation voltage, fOOO Hz 
and f 5 V, respectively, with an Andccn-Hagerling 2700A 
capacitance bridge on a compacted powder (thickness: 
~ 0.8 mm, 0: 3.6 mm). 



B. Results and Discussion 

Figure [3] shows the measured and simulated x-ray pow- 
der diffraction patterns of the sample of CuCrC>4 used for 
all subsequent characterization. The refined atomic pa- 
rameters and the lattice parameters are summarized in 
Table |TT] and were found to be in good agreement with 
the previously published single crystal results ^ 
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FIG. 3: (Color online) (o): Measured x-ray diffraction pattern 
of CuCr04 (wavelength 0.709 A Mo-Kai radiation). Solid 
(red) line: Fitted pattern (R p = 3.42 %, reduced \ 2 = 1-15) 
using the parameters given in Table HT1 Solid (blue) line (off- 
set): Difference between measured and calculated patterns. 
The positions of the Bragg reflections used to calculate the 
pattern are marked by the (green) vertical bars in the lower 
part of the figure. 

Figure @] displays the optical spectrum of CuCr04 
which is consistent with the deep brownish-red color of 
the CuCr04 powder. The spectrum is dominated by 
a strong absorption band centered at 2f500 cm -1 (466 
nm) which we attribute to an 2 ~ — > Cr 6+ charge trans- 
fer transition, in agreement with observations for other 
hexavalent chromates . 40 ' 41 In the near infrared regime 
(NIR) the spectrum exhibits a maximum at = f 3000 
cm -1 with a tail extending down to <~ 7000 cm . Two 
subsequent faint shoulders are seen within the slope at 
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v 2 = HOOO cm" 1 and v x ~ 8000 cm" 1 . 




FIG. 4: (Color online) Powder reflectance spectrum of 
CuCrOt. Black vertical bars mark the ligand-field transi- 
tion energies, for the CuOg distorted octahedron, obtained 
from AOM calculations. We show the Kubelka-Munk rela- 
tion, (l-R f ) 2 /(2 R f ), where R f = J(CuCr04)//(BaS0 4 ) and 
/(CuCrCU) and /(BaSO-i) are the reflected light intensities of 
the sample and the BaSCU standard, respectively.— 

Using ligand-field considerations (see below) the ob- 
served absorption bands v 2 and £3) can be assigned 
to Cu 2+ d-d transitions, 2 B\ g — > 2 A lg (z 2 — > x 2 - 
y 2 ), 2 B lg 2 B 2g (xy -> x 2 - y 2 ), and 2 B lg 2 E g 
(xz,yz — > x 2 - y 2 ), respectively^ 3 - - — 

From &i, i>2, and £3 the crystal field splitting, lODq for 
CuOg, can be calculated using the relation, 

10^ = *>3-(^p)-(|), 

which yields a value of lODq ~ 8300 cm -1 . This value 
is similar to crystal field splitting values previously re- 
ported e.g. for Cu 2+ aquo-complexes4£ 

UV/vis spectra for CuCrC^ have been reported before 
by Baran and an assignment of the observed transitions 
has been been proposed^ Based on calculations within 
the framework of the angular overlap model (AOM)^r— 
we argue that this assignment has to be revised. 

Within the AOM model the pairwise interactions of 
the ligands with the d-orbitals are encoded into the pa- 
rameters, e<7, en, x and e Vty which take care of interac- 
tion along and perpendicular to the Cu - 0» (i = 1, . . . , 
6) bond, respectively. The energies of the individual d- 
orbitals are obtained by summation over all pairwise in- 
teractions. The variation of the AOM parameters e a i 
with the Cu - Oi distance has been taken care of by, 

e al oc 1/r™. 



An exponent of n w 5 is derived from electrostatic and 
covalent theoretical bonding considerations! 42 ' 48 ' 49 Mea- 
surements of the pressure dependence of IQDq pointed to 
a similar exponent 5 < n < 6^ For the sake of simplicity 
we have chosen e^ %x = e„ t y = l/4e CT . AOM calculations 
have been performed using the program CAMMAG i 51 ' 52 
Table IIIII summarizes the parameters which have been 
used for these calculations. The resulting transition en- 
ergies marked by vertical bars in Fig. [4] are in good agree- 
ment with the centers of the experimentally observed ab- 
sorption features. 



Ocq Pax 

d (Cu-O) (A) 1.965 (4x) 2.400 (2x) 
e a (cm -1 ) 5600 2061 

e„,x (cm -1 ) 1400 515 

e T|B (cm -1 ) 1400 515 



TABLE III: Parameter used in the AOM calculations. The 
equatorial plane forms a rectangle with the equatorial O eq - 
Cu - C q bonds enclosing an Z of 81.92° and 98.08°, respec- 
tively. The Racah parameters amounted to B — 992 cm -1 , 
C — 3770 cm -1 yielding a ratio C/B — 3.8, as given for the 
free Cu 2+ ion.— As for the aquo-complex the nephelauxetic 
ratio j3 was chosen to be 0.80, and the spin-orbit coupling 
parameter £ = 664 cm" 1 was reduced by 20% as compared 
to the free ion value . 43 ' 53 

In addition to the energy of the excited electronic 
states of the isolated CuOg unit, its magnetic proper- 
ties are also obtained from the AOM calculations. The 
paramctrization leads to an average g av = 2.18 and a 
strongly anisotropic g-tensor with g x = 2.07, g y = 2.07, 
and g z = 2.39 along the principle axes. The z-direction 
of the g-tensor lies along the Cu - O ax bond direction. 

The results of the specific heat measurements for mag- 
netic fields of T and 9 T are displayed in Fig[5] The T 
data reveal a rather broad, smeared, A-type anomaly cen- 
tered at 8.2(5) K marking the onset of long-range mag- 
netic ordering. Within experimental error the data mea- 
sured in a magnetic field of 9 T are identical to those 
obtained at T. The plot of C p /T versus T given in the 
low right inset of Fig. [5] enables the estimation of the 
entropy contained within the anomaly, which equates to 
~ 0.6 J/molK or ~ 10 % of the expected entropy of a 
5 = 1/2 system, i?ln(2), where R is the molar gas con- 
stant. 90% of the entropy has already been removed by 
short-range afm ordering above Tn- 

At low temperatures, the heat capacity comprises of 
a phonon and magnon contribution. The temperature 
dependence of the phonon contributions to the heat ca- 
pacity can be described by a Debye-T 3 power law. The 
magnon heat capacity at low temperatures varies with a 
power law depending on the spin wave dispersion rela- 
tion and the dimensionality of the lattice. For a three- 
dimensional (3D) magnetic lattice, one obtains a T 3 
power law for afm magnons, and a T 3 / 2 power law for 
fm magnons M. The C p /T 3 /' 2 versus T 3 / 2 plot shown in 
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the upper left inset of Fig. [5] demonstrates that at low 
temperatures the heat capacity conforms well to a T 3 / 2 
power law, with the coefficient of the fm magnon contri- 
bution given by the non-zero intercept with the ordinate, 
7, according to, 



Cp/T 3 / 2 = /3T 3/2 + 7 



3/2 



(2) 



where /? is related to the Debye temperature, #d (0) at 
zero temperature via, 



5 6> D 



(3) 



with M = 6 being the number of atoms per formula 
unit of CuCr0 4 . While 7 can be expressed as^ 



B_n3/2 



(4) 



Here we have assumed that the Cu ribbon chains are 
coupled to neighboring chains by J±, which we associate 
with the fm inter-chain spin exchange constant J2 (see 
Table HI. 

By using eq. ([3]) we ascertain 9d(T — > 0) to be, 



9(T -> 0) = 138(3) K, 
and from the intercept and eq. ([2]) we obtain 7 as 

7 = 1.03(2) x 1CT 2 J/molK 5/2 . 

By using J± ~ J2 ~ 12 K (see below) and 5=1/2, the 
pre-factor A in eq. ^ amounts to 

A w 0.15 J/molK. 

Figure [6] summarizes the results of our EPR measure- 
ments. Near 3.4 kOe a single rather broad (peak-to-peak 
linewidth AH pp w 0.8 - 1 kOe) symmetric resonance line 
was observed. It can be well fitted to the derivative of a 
single Lorcntzian absorption line with a small contribu- 
tion I a I < 0.04 of dispersion according to 



dR 



abs 



d AH + a(H-H IC 



AH + a{H + H res ) 



dH 



dH {H 



H ) 2 



AH 2 (H 



H ICS ) 2 + AH 2 ■ 
(5) 

As the linewidth (half-width at half-maximum 
(HWHM)), AH, is of the same order of magnitude as 
the resonance field, H ICB (see Fig. H{a)), in eq. ([5]) we 
took into account both circular components of the excit- 
ing linearly polarized microwave field and therefore also 
included the resonance at negative magnetic fields cen- 
tered at —H ICS . 




FIG. 5: (Color online) (Black) o and (red) A: Heat capacity 
of CuCr04 at T and 9 T, respectively. The latter data have 
been shifted by +0.5 J/molK. Upper left inset: C p /T 3/2 plot- 
ted versus T 3 ^ 2 to highlight the low-temperature T 3 / 2 power 
law. The (red) solid line is a fit of the data to eq. ((2j with pa- 
rameters given in the text. Lower right inset: C p /T depicted 
against T in the low-temperature regime. 



The resonance field of the room temperature powder 
spectrum corresponds to a (7-factor of 2.117(2). Upon 
cooling a slight increase of the ^-factor with saturation 
to a value of ~ 2.125 below 150 K was observed (Fig. 
|6fc)). Such a value is somewhat lower than the expected 
average value g av ascertained from the AOM calculations. 
The resonance line is too broad to resolve the anisotropic 
g-factors which range between ~ 2.39 and ~ 2.07 (see 
above). 

The integrated intensity of the EPR resonance, I(T) 
which is proportional to the spin-susceptibility, increases 
with decreasing temperature down to ~ 15 K where a 
hump occurs. Above ~ 150 K, I(T) follows a Curie- Weiss 
type temperature-dependence, 



J(T)cx 



T — 



EPR 



(6) 



with 



9epr « -60(5) K. 



The negative T-axis intercept indicates predominant 
afm spin exchange interactions. Deviations from the 
Curie- Weiss type temperature-dependence are ascribed 
to short-range afm correlations, which start to develop 
below ~ 150 K, similar to the behavior of the dc magnetic 
susceptibility (see below). The decrease of the integrated 
intensity below ~ 15 K signals the onset of long-range or- 
dering. 

The magnitude and temperature-dependence of the 
EPR linewidth, AH, are similar to those observed for the 
inorganic spin-Peierls system CuGeOa or the frustrated 
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afm ID system LiCuV04^^ The linewidth exhibits a 
concave temperature dependence with a linear increase 
at low temperatures and for T — >• oo one extrapolates a 
saturation value of ~ 1.4 kOe. 

If we assume that the temperature-dependant broad- 
ening of the EPR resonance line is due to anisotropic or 
antisymmetric components in the exchange Hamiltonian, 
the constant high-temperature value can be estimated 
from the Kubo-Tomita limit as^ 



AH(T -> oo) 



1 8 2 



(7) 



where 8 indicates the deviations from the symmetric 
Heisenberg spin exchange and J is the afm symmetric 
intrachain exchange. If for CuCr04 we associate J with 
the nn spin exchange, ~ 60 K (see below), we can esti- 
mate a 8 of ~ 3 K, i.e. 5% of the symmetric exchange. 

The linear slope of the linewidth at low temperatures 
can be explained using the formulism put forth by Os- 
hikawa and Afflec k 59 ! 60 predicting 



AH(T) cx — T. 



(8) 



We find a linear slope, indicative of ID afm system, 
of ~ 2.5 Oc/K, similar to that observed for CuGcC>3 
(~ 4.5 Oc/K)^ 6 . 

The magnetic susceptibility of a polycrystalline sam- 
ple of CuCr04 was measured in magnetic fields of 1, 3, 
5 and 7 Tesla. Above ~ 20 K the susceptibilities are 
independent of the magnetic field indicating negligible 
ferromagnetic impurities. The susceptibilities, Xmoi(T), 
above ~ 150 K follow the modified Curie- Weiss law, 



Xmo\(T) 



C 



T-Q 



+ Xdia + XVV- 



(9) 



C is the Curie constant pertaining to the spin suscepti- 
bility of the Cu 2+ entities, C = N A g 2 /Lt| S(S + l)/3fc B . 
Xdia refers to the diamagnetic susceptibilities of the elec- 
trons in the closed shells, that can be estimated from the 
increments given by Selwood, which equates to -62 x 10 -6 
cm 3 /moli£i 

At high temperatures, T > 150 K, we fitted the mo- 
lar susceptibility to the aforementioned modified Curie- 
Weiss law (eq. H}. We found best agreement with the 
following parameters: 



g = 2.17(2) and 6 = -56(1)K 



and 



Xdia + Xw « +20 x 10 6 cm 3 /mol. 

This puts the Van Vlcck contribution to w +80 x 
10 _6 cm 3 /mol which is in reasonable agreement with 
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FIG. 6: (Color online) Results of the EPR measurements 
on a polycrystalline sample of CuCrC>4. (a) (o) Inverse of 
the integrated intensity. The (red) solid line is a fit of eq. 
© to the high temperature data (T > 150 K). (b) (o) The 
fitted half-width-at-half-maximum (HWHM) versus temper- 
ature, (c) (o) (/-factor versus temperature, (d) (o) EPR spec- 
trum of CuCrC>4 measured at RT with ~ 9.45 GHz versus 
applied magnetic field. The (red) solid line represents the fit- 
ted derivative of a Lorentzian absorption line (eq. [5} to the 
measured spectrum. 



what has been found for other Cu 2+ compounds (see 
Ref. [TU and refs. therein). The fitted g- factor is in good 
agreement with optical spectroscopy and the Curie- Weiss 
temperature is negative and in accordance with #epr - 

Below 150 K there are deviations from the Curie- Weiss 
law attributed to increasing afm short range correlations. 
The susceptibility passes through a broad shoulder with 
a subsequent kink at ~ 8 K whereupon it becomes field 
dependent, with a tendency to diverge for small fields. 
With increasing fields the divergence is suppressed and 
the kink becomes more apparent. By 7 T a pronounced 
rounded hump with a maximum at 14.2(2) K and a sub- 
sequent dip at 8.0(5) K become clearly visible. 

In general, GGA+J7 calculations overestimate the spin 
exchange constants typically by a factor up to 4, in our 
case 2j 34 ' 35 i 62 By taking this into account and by us- 
ing a mean field approach one calculates, from the spin 
exchange parameters summarized in Table [Til a (nega- 
tive) Curie- Weiss temperatures ranging between -38 K 
to -45 K, consistent with the experimental observations. 
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Our GGA+U calculations indicate that CuCr04 can 
be described by a Heisenberg ID chain with afm nn 
and afm nnn spin exchanges, with significantly weak 
inter-chain interactions {J-z/Jim < 0.19). Therefore, 
we modeled the magnetic susceptibility of CuCr04 
against exact diagonalization results for the susceptibility 
Xchain(ff, J n na) of a single chain provided by Heidrich- 
Mcissner et a^,, 63 i 64 with 



Jnn/ Jn 



(10) 



Interchain spin exchange is treated within a mean- field 
approach according to,— 




100 200 300 

r(K) 



7(K) 



Xmol(T) 



Xchain(^) 

1 - A 

Xchain 

(T) 



+ Xo- 



(11) 



By using the already known values, %o = Xdia + Xvv 
= +20xl0~ 6 cm 3 /mol as found from the fit of the high 
temperature magnetic susceptibility and a g-factor of 
2.13 obtained from the EPR measurements, the simu- 
lated results can be compared to experimental data. The 
mean-field parameter, A, in cq. (|llj) can be ascribed to 
the inter-chain spin exchange interactions according to^ 



A = (zi Ji + z 2 J 2 )/N A g 2 ^, 



(12) 



wherein, z± = 4 and z 2 = 2 count the number of spin 
moments with which a chain spin interacts through the 
inter-chain spin exchange interactions, J\ and J 2 , respec- 
tively. Guided by the GGA+U results, the ratio a is pos- 
itive and in the regime of 1.5 to 2.5. Within this range 
for a we find best agreement of our experimental data 
with the model calculations for, 



a « 2, implying J mm = -27(2) K, 
and a positive A, which amounts to 



FIG. 7: (Color online) (main panel) (o) Temperature depen- 
dence of the molar magnetic susceptibility, Xm, taken at 7 T. 
Colored solid lines represent the exact diagonalization results 
by Heidrich-Meissner et al. for various ratios of J nn /Jnnn, 
1.5, 1.75, 2 (red solid line), 2.25 and 2.5, from top to bottom, 
respectively. See text for more details. The dashed line is 
the magnetic susceptibility of a 5=1/2 Heisenberg chain with 
afm uniform nn spin exchange of -27 K.— (a) red symbols: 
heating data, blue symbols: cooling data. Xmoi versus tem- 
perature for various magnetic field, (b) (o) Reciprocal molar 
susceptibility versus temperature with a fit ((red) solid line) 
to a modified Curie- Weiss law (eq. ((§])). 



The inter-chain spin exchange can also be estimated 
from the Neel temperature, Tn, which, according to the 
heat capacity data, amounts to (see above); 



.2(5)K. 



Yasuda et al. calculated the Neel temperature of a 
quasi ID Heisenberg antiferromagnet on a cubic lat- 
tice with the isotropic inter-chain coupling Jj_, inducing 
3D long-range magnetic ordering at a Neel temperature, 



A = 7(l)mol/cm 3 . 

Figure [7] shows a comparison of the measured data and 
the mean-field corrected exact diagonalization results. 

A > indicates that the dominant inter-chain spin 
exchange is fm, consistent with our density functional 
calculations. The DFT calculations indicate Ji w - 
l/4xj 2 , irrespective of U e g. From eq. (jT2"j) using 
A = 7(1) mol/cm 3 we derive a value for J 2 which amounts 
to 



J 2 = 12(2)K. 

This value is in good agreement with the scaled DFT 
result, sec Table fl] 



T N /|J±| = 0.932 \jln{A) + Unln{A), (13) 

where A— 2.6 J\\/T-^ and Jy is the intrachain spin ex- 
change constant. If we assume J\\ to be our J nn ~ -60 
K we find the inter-chain coupling to be 

|Jx|«5K, 

consistent with the value obtained, from A. The dif- 
ferences may arise, since our real system has two differ- 
ent inter-chain coupling constants, J\ and J 2 , as indi- 
cated by our density functional calculations. Addition- 
ally, CuCr04 has a nnn spin exchange J nnn , which is not 
included in Yasuda's model. 
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7(K) 



FIG. 8: (Color online) (a) Colored symbols represent the rel- 
ative dielectric constants, e r , versus temperature for different 
applied magnetic fields, as given in the legend, (b) The zero 
field relative dielectric constant is shown by the solid black 
line within a greater temperature range, (c) (o) The relative 
dielectric constant versus applied magnetic field at a temper- 
ature of 5.2(1) K. 

Figure [5] displays the temperature and magnetic field 
dependence of the relative dielectric constant, e r , of a 
compacted polycrystallinc sample of CuCr04. 

At room temperature, a value of ~ 48 was found for e r . 
With decreasing temperature, e r is seen to decrease in a 
smooth fashion, until it passes through a shallow double 
maximum between 35 and 15 K, possibly indicating some 
magnetostriction induced by short range magnetic order- 
ing processes above Tnt (see Fig. [8] inset (b)). At 10 K a 
value of e ~ 4.35 was measured. Long-range magnetic 
ordering leads to a sizeable ME effect as evidenced in 
the e r , however, with a rather broad anomaly extending 
over the whole temperature range down to 3 K. Indica- 
tion for a sharp spike near Tn, as is frequently found in 



multifcrroic systems, has not been seen. Similar broad 
anomalies, originating at Tn, have been in seen in CuCl2 
and CuBr2 - 68 ' 69 In zero field a steep increase of e r is seen 
to occur below ~ 8.5 K with a broad slightly asymmet- 
ric hump centered at ~ 5.35 K. In zero field the increase 
of e r from the paramagnetic phase to the maximum of 
the hump amounts to ~ 6%. Applying a magnetic field 
decreases the ME anomaly and moves the maximum to 
higher temperatures. The onset of the ME anomaly is 
not seen to move, in accordance with the aforementioned 
C p measurements (see Fig. [5J inset (c)). The decrease of 
e r with a magnetic field starts above ~ 1 T and tends to 
saturation at sufficiently high fields. 

In summary, CuCrCU represents a new ID quan- 
tum antifcrromagnet with a remarkable pronounced ME 
anomaly below the Neel temperature of 8.2 K. Our den- 
sity functional calculations indicate that, to a first ap- 
proximation, the spin lattice of CuCr04 is a ID Heisen- 
berg chain with the unique situation that both, nn and 
nnn, spin exchanges are afm. J nn />/nnn is found to be 
close to 2, which places CuCrC>4 in the vicinity of the 
Majumdar-Ghosh point. The presence of sizeable fer- 
romagnetic inter-chain spin exchange interaction leads 
to long-range magnetic ordering. The occurrence of the 
rather large ME anomaly below the Neel temperature is 
taken as evidence for non-collinear, possibly helicoidal, 
spin ordering in the ID chains. CuCr04 therefore rep- 
resents a new interesting example for an unusual type-II 
multifcrroicity system. Neutron scattering investigations 
are scheduled to clarify the exact nature of the magnetic 
ground state of CuCrCU. 
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